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Abstract 

The CDF experiment has reported an excess of high-p^ jets compared to previous 
next-to-leading order QCD expectations. Before attributing this to new physics effects, 
we investigate whether these high-p^ jets can be explained by a modified gluon distri- 
bution inside the proton. We find enough flexibility in a global QCD analysis including 
the CDF inclusive jet data to provide a 25-35% increase in the jet cross sections at the 
highest pt of the experiment. Two possible sets of parton distributions are presented, 
and the effects of these on other existing data sets are presented. Further theoretical 
and experimental work needed to clarify unresolved issues is outlined. 
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Jet production in hadron collisions at the Fermilab Tevatron is an important process 
which presently provides the highest energy for studying hard scattering dynamics. Quark 
substructure or other new short distance physics would, if present, change the cross section 
for high-pt jet production Such a deviation compared to next-to- leading order (NLO) per- 
turbative Quantum Chromodynamics (QCD) calculations, based on commonly used parton 
distributions, has been reported by the CDF experiment p[ in the range 200 < pt < 420 GeV 
from 20 pb~^ of data. These data are shown in Fig. 0, indicating a clear 40% excess at pt = 350 
GeV compared to the NLO calculation 0. The points are (Data-NLO QCD) /NLO QCD 
plotted versus the scaling variable Xt = 2pt/y/s. The theory is calculated with CTEQ3M 
parton distributions [^] and fi = pt/'i. The region of excess corresponds to 0.22 < xt < 0.45. 
As more data are analysed, this excess could become more significant, as suggested by the 
first results on the dijet mass distribution from 70 pb~^ of data 0. In order to determine 
whether this enhancement constitutes a signal for new physics, it is crucial to investigate 
possible explanations within the Standard Model. 

One well-known uncertainty concerns the dependence of perturbative calculations on the 
choices of the renormalization and factorization scales. However, for inclusive jet cross sec- 
tions, this dependence is quite small (10%) and is largely independent of pt 0. Similarly, 
changes in the strong coupling «<, resulting from variations of Aqcd (such as a comparison 
between calculations using CTEQ2M and CTEQ2ML [Q, or the recent MRS0 parton dis- 
tributions) mainly affect the normalization. Another source of uncertainty is the effect of 
summing large perturbative logarithms that may be important at large-x* and have been 
shown to be significant for high-mass lepton-pair production]^. A corresponding study for 
jet cross sections has not yet been carried out. In addition, the long-standing disagree- 
ment between NLO QCD and the jet Xt scaling result from CDF points to a potential 
inadequacy in the NLO calculations, or a possible mismatch between the theoretical and 
experimental jet definitions. However, it is not clear whether this effect, even if it is real, 
will extend to the Xt-region under consideration. 

Finally, there are the parton distributions which play a crucial role in determining the 
perturbative QCD "predictions" of the jet cross section. For the Xt range in question, more 
than 50% of the jet cross section is due to quark-quark scattering, and the quark distributions 
are well determined by the precise data from deep inelastic lepton-nucleon scattering (DIS). 
On the other hand, although the gluon distribution is small in this region, its contribution to 
the cross section (mainly through the gluon-quark scattering processes) is still substantial — of 
the order of 25-50%. The DIS data do not constrain the gluon much at large x, that role being 
usually played by direct photon production data in most modern global analyses.]^, ^ In 
the light of current theoretical and experimental uncertainties on direct photon production, 
it remains an open question whether the usual gluon distributions can be modified in the 
relevant x region to accommodate the observed high-p^ jets. The purpose of this paper is to 
report on a quantitative study addressed to this particular question. 

Global QCD Analysis Incorporating CDF Jet Data We have carried out a global 
QCD analysis |Q incorporating for the first time the CDF inclusive jet data. Particular 
attention is given to the compatibility of the inclusive jet data with the collection of data 
sets used in previous global analyses within the NLO QCD framework ||. 
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As seen in Fig. 0, jet production data with pt > 200 GeV {xt > 0.22), where the ex- 
cess occurs, have comparable experimental systematic and statistical errors, whereas below 
this value the systematic errors dominate. Thus, these two regions are given separate at- 
tention. For the lower pt range, although the experimental measurements extend down to 
Pt = 15 GeV, we chose to include only data with pt > 75 GeV in this study due to a number 
of potential theoretical and experimental problems relating to low-pt jets. These include: 
1) possible problems in the match between theory and experimental jet definitions, such as 
fragmentation products outside the jet cone, 2) definitions of the "underlying event" coming 
from the proton-antiproton remnants, 3) scale uncertainty of NLO QCD calculations which 
becomes non-neghgible at low pt, and 4) kt broadening (discussed later for direct photons). 
All of these affect \ow-pt jets much more than high-p^ jets. 

Our systematic study reveals that there is enough flexibility in the NLO QCD global 
analysis framework to enhance the theoretical cross section for the highest pt inclusive jets 
by 25-35% above the previous calculations. We will describe two sample parton distribution 
sets which illustrate two slightly different ways that the overall flt can be accomplished. 
The flrst, designated as the norm=1.0 jet-flt, flxes the CDF jet data normalization at the 
nominal value so that the high-pt excess points are accommodated without an overall down- 
ward shift of experimental points. However, without flxing the normalization of the jet data, 
the global analysis prefers a relative downward shift of the CDF data with respect to theory. 
The second example, the norm=0.93 jet-flt, is chosen to represent this possibility. Both 
solutions give good flts to the other data sets included in the global analysis. 

Fig. |l| includes two curves corresponding to NLO QCD calculations using parton distri- 
butions from the two new flts along with the CDF jet data: the solid line for the norm=1.0 
jet-flt, and the dashed line for the norm=0.93 jet-flt (divided by 0.93). The two new flts lie 
virtually on top of each other. The total fo^^ the 1147 DIS, Drell-Yan, direct photon, and 
CDF jet data points in the norm=1.0(0.93) jet-flt is 1160(1130), clearly quite good. Both of 
the new flts remove much of the excess of the large pt jet data, with a x^/ H^pts = 1.36, which 
is quite acceptable considering this ignores the systematic uncertainties in the jet data. The 
quadratic sum of eight different CDF systematic uncertainties is shown as a shaded band 
below the data points. While the size of the band appears independent of jet pt, the eight 
individual uncertainties are not; they must be folded in for a proper analysis of errors after 
detailed information becomes available from CDF. But for our purposes of determining if 
the jet data can be accommodated within QCD uncertainties, the proper procedure is to 
only flt the jets with statistical uncertainties and give this data set more weight in the global 
flt, then look closely at the other data sets in the flt to see if discrepancies arise. This does 
not imply that one obtains the best estimate of the true parton distributions in nature; it 
does prove that viable parton sets exist. The detailed comparison of the jet-flt partons with 
other data sets will be described later. 

The gluon distributions from the two new flts are compared with that of CTEQ3M in 
Fig. at /i = 150 GeV, which corresponds to the middle of the high-pt data range with 
fi = Pt/^. In Fig. x'^G{x) is plotted against logx. (Since x^f{x) ■ dlogx = xf{x) ■ dx is 
the momentum fraction within dx, each curve in this plot directly depicts the distribution of 
momentum fraction carried by the gluon.) In Fig. the ratio of the jet-flt gluons to that 
of CTEQ3M is shown over the x > 0.1 range. For the norm=1.0 flt, we see a signiflcantly 
increased G{x) in the large- a; region, with a compensating decrease in the medium-x region 
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and little change in the lower-a; range. For the norm=0.93 fit, G{x) is uniformly shifted down 
from the norm=1.0 fit in the range 0.05 < x, with a compensating increase in the small-x 
region. This shows that the jet data used in the fit constrain the shape of G{x) in the region 
0.08 < a; < 0.45. Also shown in Fig. |]b, is the ratio for /i = 5 GeV, which is relevant for 
discussing comparisons with direct photon data later. We also note that the as{Mz) values 
for the new fits are slightly higher, 0.116 compared to 0.112 for CTEQ3M. 



Comparisons to Deep-Inelastic Scattering and Direct Photon Data Deep-inelastic 
scattering data are indirectly sensitive to the gluon distribution through NLO corrections 
and scaling violations. But at large-x the effects on F2 due to a modified gluon distribution 
can be easily compensated by small changes in the quark distributions and in Aqcd- A 
detailed look at the shifts in F2 for the various DIS experiments shows no changes of more 
than 2% between CTEQ3 and the jet-fit results for all values of x and Q^. 

Fixed target direct photon data have usually been regarded as the main source of con- 
straint on the gluon distribution at large-x. However, the constraint is weakened if the 
theoretical uncertainties unrelated to parton distributions are significant. We will now re- 
view the relevant theoretical issues, then evaluate the effect of the jet-fit gluons in light of 
these uncertainties. For definiteness, we shall use the most widely used WA70 data as the 
point of discussion, although the results are independent of the specific experiment. 

The two most significant theoretical uncertainties for fixed target direct photons are the 
factorization scale dependence, and the possible kt broadening effect. The latter is suggested 



by a recent global study of all direct photon data |T^; it involves the likelihood that NLO 
QCD does not contain enough of the kt of the initial state gluon radiation, thus leading to 
an underestimate of low-pj photon cross sections. Since the publication of Ref. 0], there 
have been two developments which further support the basic idea of kt broadening in direct 
photon production: 1) a new calculation of collider direct photon production incorporating 
NLO QCD hard scattering plus initial state parton showers |]Tl| shows good agreement with 
the shape of the CDF direct photon data |jl2[ . 2) The preliminary, high-statistics, E706 direct 
photon data (the most precise measurement yet at fixed target energies) also shows a 
significant excess of photons compared with NLO QCD calculations. This excess is largest 
when the pt slope of the data is greatest, which again is consistent with the expectations 
from a kt broadening effect. 

In Fig. |]a, the WA70 direct photon data is compared to NLO QCD for a variety of 
scales, using conventional ABFOW parton distributions |T^. The change in theoretical 
value in going from optimized fi (used by ABFOW and MRS) to fi = pt is about 50%. This 
large variation due to scale changes provides a measure of the theoretical uncertainties due 
to higher order corrections. Next, to show the effect due to a possible kt broadening, we 
also include in Fig. Ha a curve corresponding to a scale choice of fi = pt plus an average kt 



broadening of 0.9 GeV using the algorithm of reference |]15|. The number 0.9 GeV comes 
from the WA70 analysis of their diphoton measurement [16|. We see that the broadening 
correction is also about 50%, and brings the n = Pt curve into agreement with the data. 

Even with the large theoretical uncertainties described above, one might still expect the 
fixed target data to rule out one or both of the jet-fit gluons because, naively, the differences 
between the jet-fit and conventional gluons at a typical = 150 GeV might be significantly 
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amplified at the low fi value (2 GeV to 6 GeV) of the fixed-target experiments. But this 
is not the case: the crossing point between the jet-fit and CTEQ3M gluons occurs around 
X ~ 0.4 at = 5 GeV, as shown in Fig. |^. The comparison of the jet-fit results with the 
WA70 direct photon data is shown in Fig. ^o. In the solid and dashed curves we have used 
a scale of /i = pt/2 with no kt corrections. In the dotted curve we used fi = pt and a kt 
broadening of 0.9 GeV. All three curves are consistent with the WA70 data which has a 10% 
normalization uncertainty. These results clearly demonstrate that given the uncertainties 
with scale choice and kt broadening, the new gluon distributions are fully consistent with 
the WA70 data. As mentioned above, similar results hold for other fixed target direct photon 
data sets. 



UA2 Inclusive Jet Data Of considerable interest to our study of high-Xj jets is the earlier 
UA2 inclusive jet cross section |17] measurement. The data have high statistics, are in the 
same x range as the CDF measurement, and cover a similar rapidity range. Although the 
two experiments are at different scales set by the respective pt ranges, the QCD evolution 
between the two is not significant, hence they essentially probe the same parton distributions. 
There are some important differences however. For the same Xt, the UA2 jets are at lower 
Pt, and may be subject to the additional low-p^ uncertainties that were discussed above. 
In addition, the UA2 data are based on a jet finding algorithm that less closely follows the 
infrared-safe "Snowmass" algorithm in fact the UA2 pubhcation itself expresses caution 
concerning comparisons with NLO QCD. In our NLO theory calculations for CDF we use the 
Snowmass algorithm with R = 0.7 at the parton level, while we model the UA2 algorithm 
with the modified Snowmass algorithm with R = Rsep = 1.37. Both the \ow-pt and jet 
algorithm effects warrant further study. To account for them at present, we would assign a 
larger theoretical uncertainty, 20% compared to the nominal 10%, for the UA2 calculations. 

Fig. ^ shows the CDF and UA2 jet data compared to these NLO QCD calculations 
using CTEQ3M parton distributions, and fi = ptf^. The CDF data points have statisti- 
cal uncertainties only, while the UA2 points include statistical and p^-dependent systematic 
uncertainties (this is the way the two different groups present their data). There is an ad- 
ditional 32% normalization uncertainty in the UA2 measurement, while the CDF correlated 
systematic uncertainty band is shown at the bottom of the plot. The UA2 data are system- 
atically larger than the theory (but within the normalization uncertainty), but in general 
there is no distinct shape difference as is seen in the CDF data. If one ignores experimental 
uncertainties, the two experiments disagree with each other. But clearly both experiments 
have complicated correlated systematic uncertainties that need to be understood before con- 
clusions can be drawn. 



Conclusions If the excess of high-p^ jets at CDF persists, it will be one of the most 
important challenges for QCD. Understanding what role experimental uncertainties and/or 
conventional theoretical sources play is crucial to understanding whether there is new physics 
present. This paper has considered in detail only one possible conventional theoretical source, 
parton distributions, especially the gluon distribution. Two examples were given which show 
that there is considerable room to modify the gluon distribution so that current inclusive jet 
data can be incorporated in an overall global NLO QCD fit. On the other hand, taken at 
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face value, the UA2 jet measurement disagrees with the CDF data and prefers the shape of 
conventional gluon distributions such as CTEQ3M. This disagreement may or may not be 
accounted for by experimental systematic uncertainties. 

It will very likely take many years and much work to resolve these issues. Beyond the 
obvious need for an independent, robust, method of measuring high-x gluons, we list five 
potential studies/calculations that are needed to understand the excess: 1) measuring the 
dijet angular distribution as an independent QCD test that differentiates between parton 
distributions and new interactions; 2) incorporating the correlated systematic uncertainties 
from CDF and UA2 jet measurements into the global analysis to determine if they are 
compatible; 3) confirming and understanding the jet Xt scaling result from CDF with a new 
lower energy run at Fermilab; 4) performing the large-x resummation for jet production 
as has been done for the Drell-Yan process, to see if the excess can be so explained; and 
5) performing a pt resummation of soft gluons in the direct photon calculation to reduce 
the uncertainties due to the scale dependence and kt broadening, hence sharpening the 
constraints on the gluon distribution due to these two complementary processes. Before one 
can claim there is new physics in the CDF jet excess, it is likely that all five of these future 
studies will be necessary. 

Note added in proof: A recent paper from the GMRS group |2^ has stated it is impossible 
to obtain parton distributions in agreement with the CDF jet data. Their attempt to modify 
the quark distributions to fit the jet data results in a of 20703 for 128 BCDMS data 
points. In contrast, the two jet-fits presented in this paper give rise to a of 1^3 and 175 
respectively for 168 BCDMS data points. It appears the reason that GMRS were not able 
to find satisfactory solutions like ours is that they did not allow sufficient flexibility in the 
gluon distribution shape. 
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Figure 1: The preliminary CDF jet data is compared to a NLO QCD calculation using the 
conventional CTEQ3M parton distributions (points), and the new parton distributions fit 
to the jet data (solid and dashed lines that lie on top of each other). 
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Figure 2: (a) The gluon distributions at fi = 150 GeV from the norm=1.0 and the norm=0.93 
jet-fits are compared to that of CTEQ3M: (b) the ratio of the two jet-fit gluons to CTEQ3M 
(see text). 
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Figure 3: The WA70 direct photon data is compared to NLO QCD calculations using con- 
ventional, ABFOW parton distributions in a). Different choices of scale are shown as well 
as the effect of adding additional kt broadening to the theory. In b) the WA70 direct photon 
data is compared to NLO QCD calculations using the two sets of jet-fit gluons. (see text) 
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Figure 4: The CDF and UA2 jet production measurements are compared to NLO QCD 
calculations (see text). 
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